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Abstract 

Jacobi is one of the most famous mathematicians of his century. His 
name is attached to many results in various fields of mathematics and 
his complete works in seven volumes have been available since the end 
of the xix"" century and are very often quoted in many papers. It is 
then surprising that some of his results may have fallen into oblivion, at 
least in part. We will try to describe some of Jacobi's results on ordinary 
differential equations and the available, published or unpublished material 
he left. We will then expose the selective interests of his followers and their 
own contributions. 

There are in fact many interrelated results: a bound on the order 
of a differential system, a necessary and sufficient condition, given by a 
determinant, for the bound to be reached, an algorithm to compute the 
bound in polynomial time, and processes for computing normal forms 
using as few derivatives as possible. 

We give for all of them the form under which they could have been 
proved or rediscovered, sometimes independently of Jacobi's findings. In 
conclusion, we give the state of the art and suggest some possible ap- 
plications of Jacobi's bound to improve some algorithms in differential 
algebra. 



1 Introduction 

In two posthumous articles |17[|19) . which have been recently translated [181 . 
Jacobi has introduced a bound on the order of a system of n ordinary differential 
equations in n unknowns. Let A := {ai,j) be the matrix such that Oij- is the order 
of the equation Ui in the unknown function Xj. Let J — maXo-es„ ^11=1 0-1,(7(1)- 
A sum X]r=i Oi^a(i) is called a transversal sum and J is the maximal transversal 
sum . He claims that: 

Jacobi's bound. — The order of the system is hounded by J. 
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The bound is still conjectural in the general case. 

Jacobi's algorithm. — Jacobi gave an algorithm to compute the bound 
in polynomial time, viz 0{n^) operations, instead of trying the nl permutations. 

It has been forgotten and rediscovered by Kuhn in 1955 flD], using Egervary's 
results (see Schrijver's paper ^54^ for historical details). The idea is to find a 
canon, i.e. A G N" such that, in the matrix (a^j- + Xi) one can select maximal 
entries in each column that are located in different rows. Jacobi's algorithm 
computes the unique canon with a minimal n-uple of integres A; that we denote 
by £. Let A = max^^^, ai = A — ii and /3j = max^ j- — a^. The truncated 

jacobian matrix V is the matrix (j)ui/ dx^"'^'^^''^ . 

The truncated determinant condition. — Jacobi claims that the or- 
der of the system is equal to the bound J iff \y\ 7^ 0. 

This implicitly assumes the strong bound, defined with the convention ord^; . = 
—00 if Mi is free of Xj and its derivatives, but he gives no detail about what should 
be done in such a case. 

The shortest reduction method. — Jacobi also asserts that it is pos- 
sible to compute a normal form using only £i derivatives of equation Ui and that 
it is impossible to compute one using a smaller number of derivatives. 

This implicitly assumes |V| 7^ 0; if not, a greater number of derivatives may be 
required. This is only generically true: for some particular systems, it is possible 
to differentiate Ui at most — ii^i times, for i < n, assuming Ai > • • • > 
A„). Jacobi also gives a bound on the order of derivation of the Ui required to 
compute a resolvent representation, using Xj as a differential primitive element — 
assuming it is one. Then, Ui must be differentiated a number of times equal to 
the maximal transversal sum of the matrix obtained by suppressing the row i 
and the column j in A. 

The aim of this paper is to describe the content of Jacobi's two papers 
[T7lll9j . the genesis of these results, the history of research on the subject and the 
state of the art. We also describe some related documents from Jacobis Nachlafi, 
kept in the Archiv der Berlin-Brandenburgische Akademie der Wissenschaften 
and give a complet€0 list of references. 

2 Unpublished manuscripts 

Jacobi himself is possibly the first to have forgotten his own work. According to 
Koenigsberger [32j . his manuscripts on this subject were written around 1836 
and were intended to be a part of a forsaken project of a great work on dif- 
ferential equations. Part of it was incorporated in his long paper on the last 
multiplier [21j, but the bound itself was never published in his lifetime. The 
many versions of the text, containing numerous corrections, suggest that Jacobi 
was not satisfied of the redaction. However, these manuscripts were clearly in- 
tended for publication at the time he wrote them, as it is suggested by the many 
typographical precisions in German, written in Kurrentschrift in the margins. 

^For the best of our knowledge. . . any information about material or sources not men- 
tionned here is welcome. 
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In his preface to [VD], were the second article [12] was first published, Clebsch 
asserts that it is a little posterior to Jacobi's lectures at Konigsberg university 
during academic year 1842-43, when Borchardt was his student there. 

On page 2231 of manuscript 11/23 a), one finds a reference to a paper |22) 
published in 1834 and his work on normal form computation is mentionned in 
the second part of the last multiplier paper [^T], published in 1845, but one 
may think that the results were obtained before the redaction of that paper; 
the first page of the manuscript 1/58 a) shows that Jacobi began to correct it 
in Roma, in december 1843. Assuming that manuscripts 11/13 b), 11/22, 11/23 
a), 11/23 b) where written at the same time, as many similarities of style and 
content suggest, it must have been between 1834 and the end of 1843, certainly 
not later than 1845. 

These results are a by-product of his work on the isoperimetric problem: 
"Let U be a given function of the independent variable t, the dependent ones x, 
y, z etc. and their derivatives x' , x" , etc., y' , y" , etc., z' , z" , etc. etc. If we 
propose the problem of determining the functions x, y, z in such a way that the 
integral 

' Udt 



I' 



be maximal or minimal or more generally that the differential of this integral 
vanishes, it is known that the solution of the problem depends on the integration 
of the system of differential equations: 



" dx dt ^ dt^ 

, au ,2 dU 

arr 7 <^ 77 

= — ^ H ^ - etc 

" dy dt ^ dt^ 



— ^ §^ H n?| etc. etc. 

az dt dt-^ 



I will call these in the following isoperimetric differential equations ..." (see 
Jacobi's last multiplier article [H], GW IV p. 495). 

If the highest order derivative of Xi in U is x\'''\^ the order of Xj in the 
i''^ isoperimetric equation is Ci + Cj. Then, if the are not all equal to the 
maximal order e :— max^ Ci, we cannot compute a normal form without using 
auxiliary equations obtained by differentiating the i**^ isoperimetric equation Ai 
times with Aj = e — e^. It is also clear that J = 2 is equal to the order of 

the system, provided that the Hessian matrix [d^U / dxidxj) has full rank. We 
understand how this example can have inspired the whole theory. 

In a letter to his brother Moritz, on Sept. 17'^ 1836 Jacobi writes about 
a huge manuscript on mechanics: "/ came accross some very abstract ideas about 
the treatment of differential equations that appear in mechanical problems, for 
these differential equations, with their special form, allow some simplifications 
for the integration, that had not yet been remarked. These considerations will be 
all the more important, I think, as they extend to the differential equations that 
appear both in the isoperimetrical problem and the integration of partial differen- 
tial equations of the first order. " It seems that the remaining manuscripts come 
from the time of these first investigations, so between 1836 and 1840, rather 
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than from some later attempt. The tables A-H given in are written on 
page 2250. a of the manuscript; on the back there is a table of the doubles of 
prime numbers equal to 1 modulo 8, from 2018 to 20018, a material that could 
reflect the strong interest of Jacobi in number theory and prime numbers, a 
short time before the publication of Jacobi's Canon arithmeticus |24| . a table of 
discrete logarithms, in 1839. In a letter to the Academic des Sciences de Paris, 
published by Liouville in 1840 he evocates his work on mechanics. Liou- 
ville's commentaries are enthousiastic and suggest that a book will appear soon. 
But Jacobi writes to his brother in January 1841 [25J that he is embarrassed, as 
he does not have enough "breath" to achieve his huge project, that should have 
been entitled Phoronomie, the study of physical bodies motion. 

But in 1845, Jacobi had clearly still in mind to publish a study on normal 
forms computation, for he wrote: "/ will expose in another paper the various 
ways by which this operation may be done, for this question requires many out- 
standing theorems that necessitate a longer exposition." [21\ 

It is quite possible that this project was forgotten because of a change in 
Jacobi's life — who definitely left Koenigsberg to Berlin after a long trip in Italy — 
that also opened new contacts and new scientifitic issues. One may also consider 
a possible lack of practical examples for such a general method of computing 
normal forms. The algorithm may have suffered the same absence of contem- 
porary applications. Jacobi was right claiming that the problem of computing 
the bound was of interest by itself, but the economical questions that strongly 
motivated the mathematicians of the last century were not yet considered in the 
middle of the Xix'^ century. (See section [T]) 

3 The publication of the manuscripts 

Jacobi's widow gave the manuscripts he left to Dirichlet who began to work for 
their publication with his friends Borchardt and Joachimsthal. Very few docu- 
ments remain from their work and the best source seems to be Koenigsberger 
|32| . The papers were in great disorder. In order to class them, they gave a 
number to each page. These numbers appear on the envelops were pages that 
seemed to form a single document or to be related were stored. 

Borchardt entrusted Sigismund Cohr0, who worked on the publication of 
some others manuscripts of Jacobi, with the documents related to the bound. 
Cohn indentified (see 11/13 a) two sets of manuscripts suitable for publication 
11/13 b), 11/23 b) and worked on a transcription 11/13 c) of these sometimes 
hardly readable texts. After his death in 1861 [23J, the work was continued 
by Borchardt who published the first paper [17] in his journal in 1865. The 
second [19J was published by Clebsch in the volume Vorlesungen iiber Dynamik 
[VD] in 1866. This one was quoted by Sofya Kovalevskaya in one of her most 
famous articles [39 in 1875. The fact that these papers were written in latin 
did not seem to have been a trouble at that time. Cohn and Borchardt could 
easilly write themselves some paragraphs to fill gaps in the manuscripts and 

■^Almost nothing is known about him. There was a student of that name at Konisberg 
university in 1842-43; in 1846 some Sigismundus Cohn defended in Breslau an inaugural 
dissertation entitled De medicina talmudica. According to the vita that follows this work in 
the copy kept at the library of Alliance Israelite Universelle, he was a different man with no 
interest in mathematics. 
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Borchardt even tried to rewrite full passages in order to make them clearer. In 
his biography [35], published in 1904, Koenigsberger did not translate the many 
quotations in Latin, French and Italian. 

Borchardt also wrote some kind of abstracts of the two papers, which show 
that he fully understood their content and that he did not consider what he 
published as devoid of rigor. A slightly ironical quotation of Jacobi himself 
"Tam quaestiones altioris indaginis poscuntur.^ concludes the abstract of a 
part that clearly did not satisfy Borchardt 's standards and was not kept in the 
published version (see 11/13 c). 

4 Jacobi's mathematical results 

4.1 The manuscripts 

The first paper [17J begins with the exposition of the bound and the truncated 
determinant criterion. Then, the main part of the paper is devoted to a careful 
explanation of the algorithm with complete proofs. To obtain this paper, Cohn 
put together two different texts from document 11/13 b). The first 6 pages of 
the manuscript, reproduced in Colm's transcription, were not kept by Borchardt 
(see sec. 3). They are related to the different normal forms that a given system 
could have, with a quite complete description of systems of 2 equations in 2 
differential unkowns. 

The second [l9'| begins with a fast exposition of the algorithm, without any 
proof, followed by the example of a 10 x 10 matrix. Jacobi then explains how to 
compute a normal form, using as few derivatives of each variable as possible and 
how to compute a resolvent representation for some variable a;^ , using again as 
few derivatives of each equation as possible. This paper reproduces with very 
few changes a single manuscript 11/23 b). 

Cohn considered documents 11/22, 11/23 a) as unusable because they in- 
vestigate how Jacobi's last multiplier behaves when one changes the order on 
derivatives and computes a new normal form. Borchardt and him wanted to 
avoid the multiplier theory, possibly because Jacobi decided not to include this 
material in his paper on the subject. The § 17 in manuscript 11/23 a) fos 2217- 
2220 corresponds to the same paragraph in the last multiplier paper [2ljf|, but 
the manuscript considers the general situation, whereas Jacobi retreated to the 
linear case in the published version. 

Document II/4 is not related to the bound or normal form computation, 
but rather to the last multiplier theory. It seems an interesting unpublished 
paper of Jacobi on differential equations, including results such as the linear 
independence criterion given in Ritt's Differential algebra [53ij p. 34. 

4.2 The algorithm 

A square table (a^.j) being given, Jacobi calls transversal maxima numbers being 
maximal in their column, that are located in all different rows. The idea of jacobi 

^Then these questions require further investigation. 
"See Crelles 29, Heft 3 221-225 or GW IV 403-407. 
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for computing inaxo-GS„ ^27=1 ^i,<y(i) look to what he calls a cano that is 
anxn square table where one finds a maximal set of n transversal maxima. One 
starts with a table (a^j) and computes a canon by adding to all the elements of 
row i some integer £i. At each step of Jacobi's algorithm, one tries to increase 
the number of transversal maxima. The integers £i computed by the algorithm 
are the smallest ones, meaning that there is no canon derived from the table 
ai,j such that one of these integers may be smaller. 

The algorithm may be sketched as follows. We start with the preparation 
process: we add suitable integers to the rows of the table, so that each row 
possess a maximum, i.e. an element being maximal in its column. 

The second part of the algorithm starts with the prepared table the maximum 
of the first row being the first set of transversal maxima being considered. 

A set of transversal maxima being given, we will repeatedly compute a new 
set containing one more element. We reorganise the rows and columns in the 
following way: upper rows are the rows containing the transversal maxima and 
lower rows the remaining ones; left columns are the columns containing the 
transversal maxima and right columns the remaining ones. We denote by aster- 
isks the transversal maxima and the maxima (in their columns) that are located 
in the upper rows and right columns: maxima with an asterisk are the stared 
maxima. If there is some maximal element in some lower row and right column, 
we may already add it to the set of transversal maxima. 

We say that there is a path {transitum dari) from row i to row j if some 
element of row j is equal to a stared maximum in i, or if there is a path from 
row i to row z' and from z' to j. Rows of the first class are the upper rows 
containing a stared element in a right column and all the rows to which there 
is a path from them. If there is a lower row j in the first class, there is a path 
from an upper row io with a stared maximum ao in a right column to an upper 
row ii possessing an element ai equal to the transversal maxima of io, then a 
path from row ii to row i2 possessing an element a2 equal to the transversal 
maxima of ii, then a path from row 12 to row i^. . . and at the end a path from 
an upper row ir-i to the lower row ir = j containing an element equal to 
the transversal maxima in V-i. We may then replace in the set of transversal 
maxima those located in the rows jq, . . . , V-i by ao, . . . , a^-i and get a greater 
set of transversal maxima by adding a^. There is no lower series of the first 
class iff the set of transversal maxima is maximal. 

If so, we will need to increase some rows. The rows not in the first class, 
from which there is a path to a lower row form the third class. The second class 
contains the remaining rows. We will increase all the elements of the third class 
rows by the minimal integer such that one of them become equal to a stared 
element /? in a row of the first or second class. The computations of the maxima 
in each column and of the partition of the rows in classes may be done in 0{n^) 
operations. 

If the element f3 belongs to the second class, its rows goes to the third, as 
well as all the rows from which there is a path to it, so that the change in the 
classes partition may be computed in 0{n). This may happen at most 0{n) 

^The choice of this strange word may be related to its use in the title: Canon arithmeticus, 
that according to Schumacher came from a play on words: the computations were done by a 
Kanonier Unteroffizier. See Jacobi's Briefwechsel |25l note 2 p. 62. 
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times before there is no more elements in the second class. If j3 belongs to 
the first class, some lower series of the third class will go to the first, and the 
number of transversal maxima will be increased, as we have seen. So, we need 
to perform the partition in classes and exhaust the second class, with total cost 
0{'n?), at most n times before we get a maximal set of n transversal maxima. 
The complexity of the whole algorithm is 0{'n?). 

Basically, this algorithm is the same as Kuhn's "Hungarian method" (see sec- 
tion [7|), but Kuhn also adds constants to the columns, and not only to the rows. 
The complexity of Jacobi's algorithm, 0{'n?), is the same as the complexity of 
the variant of Kuhn's method given by Munkres |41) . 

Jacobi proved further that the li are minimal. This result is a consequence 
of the two following propositions: i) there is no unchanged row of the third 
class, i.e. with £i = 0; ii) the numbers added to the third class series in the 
algorithm are the minimal ones that may change the partitions in classes. 

On may remark that this is the only place in the two posthumous papers 
[T71 [TOpI where Jacobi provides complete proofs of his results. 

4.3 The shortest reduction in normal form 

Jacobi provides |19| a method to compute a normal form, using as few derivatives 
as possible of the given equations 



His results are only generically true; it is easy, as we will see, to make the 
requested hypotheses explicit. First, he implicitly assumes that the truncated 
determinant does not vanish. In order to provide rigorous statements, we need 
to translate his findings within the framework of some formalism that did not 
exist at this time. We will use here an elementary approach that relies on diffiety 
theory [49J. 

Definition 1 . — For short, we say that Xi is a canon for the order matrix 
(ciij), if {ai,j+Xi) is a canon, i.e. a table possessing a maximal set of transversal 
maxima. 

Let £i be the minimal canon of the order matrix (aij'j^, with Oij = oidxjUi, 
A = niaxili, Ui = A — £i and j3j = maxiai^j — ai. The truncated jacobian 

matrix Vu is the matrix (^dui/ dx'^"^^^'^^ . 

We will say that an ordering < on derivatives is a Jacobi ordering for the 
system [7] if ki — pi < k2 — P2 implies 2:^^^^ < a^^^^'' • 

Jacobi calls an explicit canonical form a .system 



where the functions fi only depend on derivatives of Xj smaller than Cj . 

Jacobi's shortest reduction method may be expressed in the following way. 

® Together with a short passage in the second paper (191 related to resolvent computation 
and using the same kind of combinatorial arguments. See subsection 14.41 
'^That we may compute using Jacobi's algorithm, see subsection 14.21 



. . . , 



Xn) = 0,...,u„(xi,...,a;„) = 0. 



(1) 




fi{x), 1 <i<n 
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Theorem 2. — i) Assume there exist some functions X : 1 S:j{t), [a, b] ^ 
R that form a solution of the system u, such that |V„| 7^ 0. Then, there exists 
a normal form of u for a Jacohi ordering, of which x is solution, that may he 
computed using equation Ui — and its derivatives up to the order £i . 

ii) Assume that and all its minors of order n — \ that do not contain row 
fi have full rank, then there is no normal form of u that may be computed using 
only derivatives of up to an order strictly less than . 

Proof. — i) We may reorder the unknowns Xj and the equations Ui, so that 
the sequence ai of definition [1] is increasing and the n principal minors of Vu 
have non vanishing determinants. Let be the determinant of the fc"^ minor. 
Consider the jacobian matrix of the system <i<n, 0<fc< l-i}, 

with respect to the set of derivatives E := {cc^"'''''^'''''^'' |1 <?<n, 0<fc< 
£i}. Its determinant is a product of powers of the Di and so is not vanishing. 
We can then use the imphcit functions theorem and find, on a suitable open 
set containing (a;i(0), . . . , _ _ ^ ii(0), . . . , xi^"^'^^"'^^"''), functions 

expressing the derivatives of E, depending on the derivatives of the set S := 
{xf'^ll <i<n, 0<k<at+ 13,}. Hence we get a normal form: 

of which X is a solution. The order associated to this normal form is a Jacobi 
ordering. 

ii) Assume that there exists a normal form that may be computed using 
derivatives of Ui up to order A^, and that < ifj,. Let s :— max"^j^ Ai — £i. As 
Vti ^ 0, we see that there is an element in the normal form depending of some 
it^^'"'' with Xig ~ £ig = s, and of the form Xj^'"'^'^^^^ = g{x), where A = max^^^. 
Remaining elements in the normal form have left side derivatives Xj^' +^+'* ) ^ 

with s' < s, so g does not depend on derivatives Xj^'^^^'^'^: this implies that 
the minor of Vu obtained by supressing row «o and column jq is not of full rank, 
a contradiction. ■ 

Jacobi gave no proof for these statements; the style of the article [111 and 
of most of his manuscripts on the subject is that of a mathematical cook-book: 
in the best case, proofs are reduced to a short sketch. Ritt surmised [52J that 
the bound was suggested to Jacobi by such considerations on normal form com- 
putation. This natural assumption is confirmed by Jacobi's claim [17J that his 
normal form reduction method provides an alternative proof of the bound. We 
easily see, according to the shape of the normal form above, that the order of u 

is a% + fi%^ J- 

Jacobi also claims that there are as many normal forms of this kind as 
there are permutations a such that X]r=i = J ■ This is perhaps his single 

claim that does not stand, even under suitable genericity hypotheses, as shown 
by the example x'l + X2 + x'^ = 0, x'2 = 0, X2 + X3 — 0: we only have one 
possible permutation, but two possible normal forms for shortest reductions: 
x'l = —X2 — x':^, x'2 = 0, = —2:2 and x'l — — — x'^, a^3=0, X2 = —x^. 

If we want to compute a normal form using a different kind of orderings, 
we may need to differentiate the defining equations a greater number of times. 
Jacobi provided bounds for the computation of resolvents. 
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4.4 Resolvent computation 

In § 4 of the second article [19J, Jacobi investigates the computation of resol- 
vent^ i.e. of normal forms such that one equation of order J, a;|^^ = fidxig), 
only depends on a single indeterminate and the remaining ones, Xi = fi(xig), 
i ^ iQ, express the other variables as functions of Xi^ and its derivatives of or- 
der lower than J. He writes: "As mathematicians use to consider such kind of 
normal forms before others, I will indicate how many times each of the proposed 
differential equations ui — 0, U2 — 0, . . . u„ = are to be differentiated in or- 
der to make appear auxiliary equations necessary for that reduction. " Jacobi's 
presentation requires the implicit hypothesis |V„| ^ 0, but his results stand in 
the more general situation of quasi- regular svstems|50|. 

Let A := (uij) be defined as in the introduction, Jacobi starts with the new 
matrix A' defined by adding to each row of A the corresponding numbers £i 
defined above. He assumes that a set of transversal maxima is chosen in the 
canon A'; he denotes transversal maxima with asterisks and calls them stared 
terms, and he underlines the terms being equal to the stared maximum located 
in their column. He says that row ii is attached (annexa) to row io if it contains 
a stared term equal to some underlined term in row io or in some row attached 
to row iQ. The row io is implicitely assumed to be attached to itself. If not 
all rows are attached to row io, he increases all the rows attached to it of the 
same minimal number that makes new underlined elements to appear and new 
rows to be attached to row ?o- The process is to be continued until all rows 
are attached to it. We get a new matrix A"; the last step is to increase all its 
elements of a same number that makes the stared term in row iq become equal 
to the order J of the system u: we get a new matrix A'" , which is obtained by 
increasing the rows of A by a sequence of numbers hi, I < i < n, that will be 
the researched orders of derivation. 

Example 3. — We consider a simple system x'( — x'2 — 0, X2 — x^ = and 
X3 — X2 = 0. For this system, the order matrix 



2* 


1 


— 00 


— 00 


2* 





—00 





1* 



is already a canon. It is not possible to construct a resolvent representation 
using X2 or 2:3, but it is possible with xi. There is no row attached to the first 
one. Increasing it by 1, row 2 becomes attached to row 1. Then, increasing rows 
1 and 2 by 1, row 3 becomes attached to row 1. So, we get the matrix 

/ 4* 3 -00 \ 
A" = [ -(X) 3* 1 . 
\ -00 1* / 

The order J of the system is 5 and we need to increase all terms by 1, so that 
the stared term of row 1 be made equal to J. We get the new matrix 

/ 5* 4 -00 \ 
A'" = -00 4* 2 . 
\ -00 1 2* J 

*This word is not used by Jacobi, we borrow it from Diff. alg. [53 p. 41 
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To obtain A'" , one has increased row 1 of A by 3, row 2 by 2 and row 3 by 1, so 
one needs to differentiate the first equation 3 times, the second 2 times and the 
third 1 time to compute a resolvent representation for 

Theorem 4. — Assume that | V„| ^ and that Xj is a differential primitive 
element, then a resolvent representation of u using Xj^ may be computed using 
derivatives of equation ut up to order hi . 

Jacobi gave no proof for this result; we propose for the reader's convenience 
the following elementary one, that follows Jacobi's presentation. 

Proof. — Using theorem we may assume that the system admits a 
normal form of the shape xl"'"*'' = fi{x), that can be obtained using derivatives 
of Ui up to order ti at most. Let us denote by Wi the i"^ equation of this 
normal form. We may make a more precise evaluation of the order of derivation 
requested to compute Wi^. If x^°^^°~^^^'^ appears in Uj^, then we will need the 
(^i — ^ioY^ derivative of Wi to compute Wi^. Then, the row i is attached to the 
row zq in A' . We may recursively prove that, if the (£i — ti^Y^ derivative of Ui 
is needed, then the row i is attached to the row iq in A' . More precisely, if we 
need to increase the row zq of A' hy s < — £ig , so that row i becomes attached 
to it, we only need to differentiate Ui up to order £i — iig — s in order to compute 

As |V„| ^ 0, the order of the system is J and for computing a resolvent, we 
need the first J derivatives of Xig; we shall differentiate equation Wig, which is 
of order Oigjg in Xig , up to order J — Oig^ig . At the beginning of the process and 
after each step of differentiation, if some x^^'^ t j ^ io appears in the right side, 
we may substitute to it the expression fj{x), and repeat the process until no 
such derivative appearj^. If the row i of A' is attached to row iq, then we may 
compute the derivative of order £ig of Wig using the derivatives of Ui of order £i 
at most, and so we may compute a resolvent using the derivatives of Ui of order 
at most J — {oig^ig +£ig) +£i = hi. If the row i in A' becomes attached to row io 
after this last row has been increased of e, then we may compute the derivative 
of Wig of order e + £ig using the derivatives of Ui up to order at most £i. So, we 
may get the derivative of order J — Oig^ig of Wig, requested to get the resolvent, 
using derivatives of Ui up to order at most J — Oig^ig — {e + £ig) + £i = hi, hence 
the result. ■ 

Remark 5. — It is easily seen that the maximal possible value for hi is 
J — Oigjg + £i — £ig = J — (Pig + A) ^ £ i . H tlic cquatlou Ui has order e^, 
the sum of the hi, that is the number of equations in the system one needs 
to solve to compute a resolvent, is maximal when is the order of Ui in all 
the variables: hi is at most X^i'^^^i^*' Y^7=i hi < [n — l)X]r=i 
other hand, assume that, after some reordering, the rows are listed by successive 
order of "attachment" to row jg = 1. We have hi — X]i>i '^i,* ^^"^ i > 1, 
^ Y.k>i°-k,k, so that Ya^i hi > I]i>i(« - 

^This means that a resolvent exists for that element, see Cluzeau and Hubert [2] for details. 
^"This idea appears in some unpublished manuscripts of Jacobi, e. g. in 11/23 a) p. 2217. a: 
"et simulac in dextra parte obvenit variabilis xyj differentiale fajj^^"^, eius e [u] substi- 
tuo valorem ffj]-" (Mathematical notations between square brackets have been changed to 
correspond with those used in our proof.) See our translation |20| p. 37. 
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Remark 6. — In more general situations, that is when [V„| — 0, one may 
need to differentiate the defining equations a greater number of times, that 
may also depend on the degree of the equations. But for quasi-regular systems 
[571 149j . Jacobi's bound for the hi still stands [50], 

Jacobi did not stop at this stage; he has also provided the following elegant 
version of his result. 

Theorem 7. — The order hi up to which one needs to differentiate equation 
Ui in order to compute a resolvent representation for u using Xi„ as a primitive 
element, is equal to the maximal transversal sum of the (n — 1) x (n — 1) matrix 
obtained by removing from A line i and column ig. 

Jacobi gives a complete proof of the equality between the order hi defined by 
his process and the maximal transversal sum of the above theorem. It relies on 
the same kind of argument as the proof of his algorithm. We refer to his paper 
[19] for detailfEl. 

We may remark that a nai've use of this theorem requires 0{n^) operations 
to compute the hi using Jacobi's or Munkres' algorithm [IT]. But Jacobi's 
process, described above, only requires 0{n^) operations, once the £i have been 
computed with the same complexity. Jacobi states first an efficient algorithmic 
version of his result, before giving an elegant — but less efficient — mathematical 
version. 

These sharp bounds provided by Jacobi for computing normal forms can be 
used in a straightforward way to improve many algorithms developped in recent 
years. His results on normal form are especially important for new methods of 
resolution for algebraic systems relying on the representation of polynomials as 
Straight-Line Programs [T2] that begin to be extended to differential systems 

[iiiBiiailD]. 

4.5 The bound 

The work of Jacobi suffers from the lack of a rigorous theory allowing precise 
definitions of the mathematical objects he considers. However, if one has in mind 
that his goal is to consider physical situations, that imply implicit hypotheses, 
his proof of the bound is not so weak as one may think at first sight. 

Proof of the bound. — The proof relies on three succesive steps. The 
first is to claim that one can reduce to linear equations. This, of course, will not 
stand for all systems, but for those expressing the laws of mechanics or other 
problems of physical interest, we can take this for granted. 

From a mathematical standpoint, the most general condition under which 
this may be done is the one given by Johnson [29], expressing that a differential 
system is in some way "regular'!^. Under such hypotheses, one may from a 
system Ui = build a new system Sui = 0, which is nothing else than what 
is described by Johnson as "Kahler differentials" [2H]- Let X : t i—^ x{t) be a 
solution of u. It is a quasi-regular solution of u at to if for all r € N. when 
substituting to any derivative x^'^'' the value x^j''\to), the jacobian matrix of the 
system u, . . . , u'-''^ has full rank n{r + 1). Then the jet of X at to is a regular 

^^See our translation ;20 p. 62—63. 

^^This was to be developped later by Kondratieva et al. 1351 1371 l49l . 
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point of the subspace V of J°°(R, R") defined by u and the order of u, that is 
also the dimension of V, is equal to the dimension of its tangent space, defined 
by Su. 

It is easily seen that, if for X, then X is quasi-regular. We say "quasi- 

regular", for some "singular solutions", according to Ritt's teminology, such as 
a; = for the equation (cc')^ — 4a; = 0, may satisfy the Johnson condition. 

Jacobi uses in fact a stronger result, but without proof, claiming that if the 
system u has order e and if its general solution depends on e arbitrary constants 
fli, . . . , Ce, then the set dx/dai is a basis of solutions of the linear system Su. I 
don't know if it was "well known" at that time. 

The next step is more surprising, for Jacobi claims, with no justification, that 
it is enough to consider a linear system with constant coefficients. There is a 
beginning of proof, striked out by Jacobi on page 2203. a of manuscript 11/13 b): 
"In exploring the order of a system, as one considers only the highest derivatives 
in the linear differential equations to which the proposed ones have been reduced, 
one may assume that the coefficients are constants. For, having differentiated 
the equations [Su — 0/ many times, in order to obtain new equations". . . This 
interrupted sentence suggests that the shortest reduction process [12] used to 
compute normal forms was the basic idea. 

In fact, provided that |V„| does not vanish, we only need to consider the 
coefficients of the leading derivative^^ when computing the order: derivatives 
of the coefficients will only affect smaller derivatives of the variables. If |V„| 
vanishes, assume that the equations Ui are sorted by increasing and that io 
is the smallest integer such that the first io rows of V„ are linearly dependent, 
satifying X^iLi ^i^i — 0, then we may in Su replace Suig by Y^^i=i^ Ci(5ui)^^'~'^'o), 
as Gig 7^ 0, this new system is equivalent to Su and the leading derivatives in 
Suig have been removed: so the new system has a strictly smaller Jacobi number 
J. This gives an easy proof by induction of the result. 

We could conclude from such considerations on normal forms of linear sys- 
tems, but Jacobi uses a different kind of argument for the last step of his proof. 
Having reduced his investigations to the case of a linear system with constant co- 
efficients, he looks for solutions of the form xj = cje^*. Substituting such an ex- 
pression in his linear equations, he gets a system of the form X)J=i W ~ 0' 
\ < i < n where deg_)^ Pi j — Oij. The number of possible values for A is the 
degree of the determinant which is at most maxo-gs,, Yl^=i '^i,<y(i) = J ^ with 
equality if |Vti| ^ 0. Jacobi does not consider the case of multiple solutions, 
etc. but there is no difficulties. 

We cannot know precisely how Jacobi could have detailed his demonstration. 
However, we have shown that we can design a complete proof, using elementary 
arguments, following the indications he left in his manuscripts, provided that 
we retreat to the safe ground of regular systems. For an account of the research 
on singular solutions during the first half of the xix**^ century, see the work of 
Houtain \27\ . 

^■^We call here a leading derivative of an equation, a derivative x^*^' such that k — f)j is 
maximal, i.e. with a maximal Jacobi order. 
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5 The second part of the XIX century 

The first publication of tfiese two papersfTTl [TO] in 1865 and 1866 and a new 
publication in 1890 in the volume V of Jacobi's complete works did not stimulate 
further research on the subject in Germany. The works of his continuators dur- 
ing the xix*^ century are very superficial. They did not seem to have considered 
the subject could deserve a real mathematical effort. 

Nanson in 1876 considers the linear case with constant coefficients, in 
2 and 3 variables. He rejects the idea of substituting Xi = qc^* in the equation, 
in order to obtain the order, claiming that one should first compute the order 
in order to be sure that a complete solution of that kind could be obtained. He 
proceeds heuristically, eliminating one variable after the other and bounding at 
each step the orders in each variable of the equations he gets. 

Jordan in 1883 [30J considered the non linear situation with 4 variables. 
He tried to eliminate X2, x^, Xi in order to compute a resolvent for xi, using 
arguments that only work in the most general situation. From heuristic con- 
siderations on the number of derivatives to eliminate, he established theorem 
[7] in four variables: one needs to differentiate Ui a number of times equal to 
hi, the maximal transversal sum of the matrix obtained by removing from the 
order matrix A row i and column 1. So, the order of the resolvent will be 
max^^j^ fli 1 + hi = J. 

Nanson, who referred to Boole ^ for systems of order 1 and Cournot 
for systems of two equations, did not quote Jacobi. Jordan — who did not quote 
Nanson — did not have a full view of Jacobi's work. He claimed that Jacobi had 
given an "indirect" proof and that he would give a "direct" one. 

The work of Chrystal in 1895 [3] was rigorous, but he only considered the 
easy linear case with constant coefficients — Jacobi's arguments only worked for 
all different eigenvalues. Ritt, who gave these references [S^ also referred to a 
paper by Sarminski {Communications of the University of Warsaw, 1902) that 
I was unable to find. 

6 Ritt's work 

Ritt, who was known to be fluent in many languages, certainly had a better 
view of Jacobi's results. However, a century after Jacobi wrote them, the style 
and spirit of mathematics did change. One expects rigorous proofs, but also 
intrinsic results, attached to geometrical objects. One thinks of varieties (or 
"components") and not of systems: more precisely, in his article [S2], published 
in 1935, a system means a component. It is remarkable that this change also 
concerned a mathematician like Ritt, who knew well some very applied style of 
mathematic^ and whose activity was dominated by the spirit of "classicism' 

One interest of Ritt in this subject was to secure results that could be ap- 
plied to components intersection and it is not the case for Jacobi's bound (see 

^^He worked performing computations in the Naval Observatory in Washington during his 
studies and helped to organise a computation group working for the US artillery during World 
War I [15]. 

^^In Ritt's obituary 1441 p. 310, E.R. Lorch writes "His media are complex function theory of 
the nineteenth century and differential equations. Much of his work could have been written 
a half century earlier. " 
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Diff. Alg. [53] p. 138-144). Having developped a theory that allows to charac- 
terize "singular components", he wonders if the bound stands for all components, 
including those that do not satisfy natural hypotheses of regularity: a difficult 
question that was certainly not considered by Jacobi. 

The less convincing part of Jacobi's "proof" is to go from time varying sys- 
tems to constant coefficients. Ritt's proof in the linear case [52J solves the 
problem. But Ritt, who only considers [551 [53] elimination orderings does not 
prove the necessary and sufficient condition for the bound to be reached, given 
by the non vanishing of the truncated Jacobian. Such a condition is, as we have 
seen in 14.51 more easily proved with an orderly ordering for an adapted order 
defined in this way: ord^^^Wi := oidx^Ui — l3j, with j3j defined as in def. [T] Ritt's 
proof relies on some simplified method for computing a characteristic set in the 
linear case. He proves the strong bound, using the convention oidx^Ui — — oo if 
Xj and its derivatives do not appear in Ui. 

His 1935 paper [52] concludes with a proof of the bound for any component 
of dimension zero of a system of two polynomial equations A and B in two 
variables, that is reproduced with a few modifications in Diff. alg. |53j p. 136- 
138. One may remark that in the 1935 article [51], a footnote precises that if 
one of the variables does not appear in A, the order of A in this variable is 0. 
In [52], a new footnote claims that we can in fact prove the strong bound. The 
requested modifications in the proof seem easy, but are not given explicitely. In 
1935 ^52j, Ritt refers to Gourin [13J for the following result: if a zero dimensional 
differential ideal I contains a zero dimensional differential ideal J^, then ordi < 
Old J. 

An important argument is not explicitly stated in the proof. Ritt reduces 
the situation to the case of two polynomials, A and E, were E depends only 
on Xi. He claims that E must effectively depend on that variable. For this, 
he needs to use the fact that a single equation in two variables cannot define a 
component of dimension 0, which is true by Diff. alg. [53] chap. Ill § 1 p. 57. 

In Diff. alg. [53] p. 139-144, Ritt also proved an important result. He con- 
siders irreducible differential polynomials in two variables, and investigates the 
order of the intersection of their general solutions. He first proves that Jacobi's 
bound stands if A and B have order not greater than unity. Then, he exhibits a 
family of polynomials of order r > 3 in xi and X2, the general solution of which 
intersects the manifold of xi in an irreducible manifold of order 2r — 3. So, the 
bound cannot stand for manifolds, but just for systems. 

7 The assignment problem 

In 1944, the R.A.F. tried to optimize the reaffectation of soldiers of disbanded 
units [54J. No practical solution could be used before the end of World War II, 
but this initiated the first research on the problem. It was then considered to 
optimize the affectation of n workers to n tasks, a.ij representing the produc- 
tivity of worker i if affected to task j. One looks for a maximum, with the 
constraint that two different workers must be given two different tasks. The 
Monge problem |47| may be considered as a first, continuous, example of such 
problems (how to transport earth from a given area to some other with the least 
amount of carriage). 
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It is not the place to give much details about the discovery of the Hungarian 
method by Harold Kuhn in 1955 [S3]. Anyway, it may be of interest to consider 
the situation from the standpoint of the transmission of mathematical results. 
Jacobi's algorithm was sleeping in papers written in a dead language, with 
titles that cannot be related to the assignment problem. It also seems that the 
mathematical community was not always of a great help for the practitioner 
who wanted to solve his optimization problem in a short time. It is amazing 
that trying n\ possibilities may have been considered as an acceptable solution, 
provided that their number is finite, in the middle of the XX*^ century (see 
Schrijver [S3] p. 8). For Jacobi, trying n! solution was not a solution at all. 
He claimed indeed to look for a solution, whereas we would rather say that we 
are looking for an efficient one. The efficiency issue was at that time — very 
strongly — implicit . 

One may also notice that rediscovering Jacobi's method took more than 10 
years, from 1944 to 1955, and that prominent mathematicians such as John von 
Neuman considered the problem. It could have been much longer if Kuhn did 
not translate from Hungarian Egervary's paper [TT] that allowed him to con- 
clude. Inspired by Konig, a pioneer of graph theory [35], Egervary considered 
the problem as a weighted variant of the maximal matching problem, but he 
did not give a polynomial time algorithm [31]. Possibly, Egervary could have 
contributed to the question himself, but it seems that the research was concen- 
trated in the eastern part of the world, mostly in the United States. It was also 
strongly motivated by economical and organizational issues and one may guess 
that it did not facilitate collaborations with eastern scientists during the cold 
war. Egervary heard of Kuhn's algorithm as late as in 1957 and went back to 
such matters with two papers on the transportation problem. Tragic circum- 
staces interrupted his research. If, after the war, the Hungarian Academy of 
Sciences provided him good working conditions that stimulated his work, he 
killed himself in 1958, persecuted by bureaucrat^. 

Richard Cohn is the first, for the best of my knowledge, to have made a link 
between Jacobi's work and the assignment problem [S], but the information did 
not spread into the optimization community before 2005. 

8 The second part of the XX*'^ century 

The second part of the xx*'^ century is dominated by the work of Richard Cohn 
and his students. 

8.1 Greenspan's bound 

Greenspan proved a different bound, in the framework of difference algebra [13] 
in 1959. It is easily translated in differential algebra: let rj = max"^j^ orda;^Ui 

(ri) 

and Tji be the greatest integer such that oidx < rj , 1 < j < Greenspan's 
bound is G := J27=i ''"i ~ rnaxf^j^ rji. It was proved by Cohn in 1980 [4] that the 
order of any zero dimensional component of an arbitrary differential system is 
bounded by G. One may remark that this result implies Jacobi's bound in two 
variables, and that it may be proved using an adapted version of Ritt's proof. 

^®K6nig commited suicide in October 1944, a few days before Budapest Jews were forced 
into the ghetto and their deportation began. 
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8.2 Lando's bound 



Barbara Lando proved in 1970 [32] the "weak bound" for order one differential 
systems. The "weak bound" means that if Xj and its derivatives do not appear in 
Ui, we use the convention ordx -Ui = 0. This result was translated in difference 
algebra |43) in 1972. Proving the strong bound for order one systems would 
imply the strong bound for any system, for the strong bound is compatible with 
the classical reduction of a system to order one equations. Lando's proof uses 
results on matrices of zeros and ones that are very close to some theorems of 
Konig and Egervary. 

B. Lando also proved that any order matrix A is the order matrix of a system 
for which the bound is reached. 

8.3 Tomasovic and PDE systems 

There was a attempt to generalize the bound to partial differential systems |55| , 
due to Tomasovic in 1976. Consider a system of n partial differential equations 
in n variables and m derivatives. Let "P be a component of of {u} and 



be the Hilbert polynomial of V. If the dimension of V is 0, then Om = 0. The 
Jacobi conjecture of Tomasovic states that Om-i < J, where J is defined as 
above, according to the order matrix of the PDE system. This conjecture had 
already been presented by his thesis adviser, Kolchin, in 1966 [331 ■ Tomasovic 
proved it for linear systems and for n < 2. A proof in the linear case was also 
given by Kondratieva et al. |36| . 

Tomasovic's results remained unpublished, due to his untimely death. His 
dissertation contains interesting material that requires further examination. 

8.4 Order and dimension 

In 1983, Cohn proved that the bound would imply the "dimension conjecture": 
every component defined by a system of r equations has differential codimension 
at most r. Tomasovic proved in the chap. 6 of his thesis [SS] that the dimension 
conjecture is equivalent to this one: // a system has a component of differential 
dimension 0, then its Jacobi number J is not —oo. 

We have seen in section [5] that Ritt's proof in two variables requires the 
dimension conjecture for r = 1. More precisely, Cohn's proof shows that proving 
the bound for a system of n equations implies the dimension conjecture for a 
system of equations in r < n variables. Furthermore, Cohn proved that the weak 
bound and the dual of "Bezout's bound' also imply the dimension conjecture. 

It is known that the intersection of two manifolds of differential codimensions 
r and s may contain components of codimension greater that r + s (see Diff. alg. 
[53j p. 133) and, as we have already seen above in sectior(5J Jacobi's bound may 
only be expected to stand for systems and not for manifolds. If the examples 

^'^It is defined as ^2^—i ordx^Ui. 
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given by Ritt in these two cases are clearly distinct, we may remark that they are 
both closely related to the structure of the singular place of the manifold. We 
still need to understand better such paradoxical behaviours and their possible 
connections. 

8.5 Kondratieva's proof 

As we have already seen in subsection l4.5l the case where Jacobi's linearization 
argument works corresponds to the regularity hypothesis defined by Johnson 
in order to prove Janet's conjecture [29] in 1978: the differentials du'^\ 1 < 
i < n, k G N, are linearly independent. Kondratieva et al. call such systems 
independent systems; they were able in 1982 |35| to prove the strong bound, 
using first linearization as in Jacobi's approach [T7], then Ritt's proof for the 
linearized system du. However, this result received little attention, the paper 
being written in Russian and difficult to find. A new proof, also valid for 
independent partial differential systems has been given in 2008 by the same 
authors [37J. 

8.6 Other works 

In 1960, Jacobi's strong bound was rediscovered independently by Volevich [55] 
for arbitrary linear systems, assuming that the truncated determinant does not 
vanish. One may also mention the works of Magnus, |45l HHI who refers to 
Chrystal and Jacobi. 

9 Beginning of the XXI*'^ century 

Hrushovski in 2004 [TB] proposed a proof for the bound in difference algebra. 
His method is completely different from those used so far in this field, but it 
does not seem possible, or at least not easy, to deduce from this result a proof 
in the differential case. 

See our article for a proof of the truncated jacobian condition in the 
framework of diffiety theory, assuming the regularity condition of Johnson. A 
generalization to underdetermined systems is also considered. Kondratieva et 
al. |37| provided a proof of Jacobi's bound for independent partial differential 
systems. 

In 2001, Pryce [5T] rediscovered Jacobi's shortest reduction method in order 
to provide a efficient method of computing power series solution of implicit 
differential algebraic systems. 

10 Conclusion 

One says that Jacobi once told to a student who wanted to read all the math- 
ematical literature before starting his research: "Where would you he if your 
father before marrying your mother had wanted to see all the girls of the world?" 
So, we may hope that he would have forgiven us for having forgotten some of 
his results. 

We see nevertheless that a closer look to the past may be fruitfuU. Contem- 
porary mathematicians who turn to computer algebra will find some common 
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spirit with the approach of these times, when a great faraihary with hand com- 
putation produced a still unformal but deep attention to efficiency. 

Some of the results presented here, beyond their intrinsic mathematical in- 
terest, could help producing improved bounds of complexity or designing new 
algorithms. E.g., the shortest reduction leads to the choice of an ordering for 
which the computation of a characteristic set may be easier. 
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Warning 

The bibliography is divided in four parts. The first, Manuscripts contains primary 
material, mostly Jacobi's manuscript, comming from Jacobis NachlaS, Archiv der 
Berlin-Brandenburgische Akademie der Wissenschaften. Documents are denoted 
by their archive index, e.g. [1/58 a)]. The second, Complete works contains the 
book were Jacobi's work were published. They are denoted by [GW ...], [VD], 
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were published. Their are denoted by [Crelle. . .] The fourth and last contains the 
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